Spacetimes containing slowly evolving horizons 

William Kavanae 

Department of Physics and Physical Oceanography, Memorial University of Newfoundland 
St. John's, Newfoundland and Labrador, A IB 3X7, Canada 

Ivan Bootl^l 

Departm,ent of Mathematics and Statistics, Memorial University of Newfoundland 
St. John's, Newfoundland and Labrador, AlC 5S7, Canada 
(Dated: February 4, 2008) 

Slowly evolving horizons are trapping horizons that are "almost" isolated horizons. This paper 
reviews their definition and discusses several spacetimes containing such structures. These include 
certain Vaidya and Tolman-Bondi solutions as well as (perturbatively) tidally distorted black holes. 
Taking into account the mass scales and orders of magnitude that arise in these calculations, we 
conjecture that slowly evolving horizons are the norm rather than the exception in astrophysical 
processes that involve stellar-scale black holes. 



I. INTRODUCTION 



Physics is most exciting far from equilibrium. Few would argue that laminar fluid flow is more interesting than 
turbulence or that a slowly cooling cup of tea is more intriguing than the rapidly boiling kettle of water that directly 
preceded it. That said, it is precisely the relative dullness of equilibrium or near-equilibrium systems that helps us to 
understand them so well. In turn, this understanding is the foundation on which much of our knowledge of nature is 
based. 

Though black holes arc intrinsically extreme objects, they too have equilibrium states. Traditionally these have 
been identified with the Kerr family of spacetimes which consists of non-evolving but rotating black holes sitting alone 
in an otherwise empty, asymptotically flat universe. Given that these are perturbatively stable and have been shown 
to be the unique stationary, rotationally symmetric and asymptotically flat vacuum black hole spacetimes, it is widely 
believed that in the absence of external interactions all black holes will ultimately settle down into a state that is 
closely approximated by one of these solutions (see, for example, Q for a discussion of these results). Much of what 
we know about black holes comes from a study of these equilibrium solutions and their perturbations. 

Of course, a key word in the above paragraph is "approximated" . Real astrophysical black holes don't each occupy 
their own asymptotically flat universe. Thus when considering a real black hole, one typically assumes that the rest of 
the universe has only a weak influence on the spacetime structure close to the black hole. Then, in that neighbourhood 
the spacetime should be perturbatively Kerr. In the last decade an alternative approach to this problem has been 
proposed. The isolated horizon programme (see |^ Q for reviews) explicitly studies the quasilocal conditions 
under which a black hole may be said to be in equilibrium with its (possibly dynamic) surroundings and then studies 
consequences of those conditions - as opposed to studying equilibrium spacetimes which also contain black holes. 

Such a characterization is important if one wishes to consider cases where a black hole may sometimes dynamically 
interact with its surroundings and at other times be quiescent. While the standard causal definitions of black holes and 
event horizons work well for stationary spacetimes, they encounter serious problems for these more general situations. 
A primary problem is that dynamic, causally defined black holes cannot even be precisely located by non-omniscient 
observers. However, there are secondary problems as well. A particularly serious one for those interested in equilibrium 
states follows from the area increase theorem. This says that an event horizon can never decrease in area and what is 
more its rate of expansion must always decrease. Thus, due to their tclcological definition, event horizons will grow in 
anticipation of future interactions even when there is no proximate cause for such an expansion. In particular, lack of 
growth cannot be used to characterize (temporary) equilibrium states (see [91 for a further discussion of these points). 

Following these arguments for quasilocal definitions. Booth and Fairhurst [6j recently proposed a characterization of 
near-equilibrium black holes : slowly evolving horizons. These "almost isolated" horizons include (most) truly isolated 
horizons as particular examples and in turn are special cases of Hayward's trapping horizons and are closely related 
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to standard apparent horizons as well as Ashtekar and Krishnan's dynamical horizons 0i 0, 13 (see for a detailed 
discussion of how these objects are related). 

In this paper, we shall consider several examples of spacctimcs containing slowly evolving horizons including certain 
Vaidya (section and Tolman-Bondi (section llV|l spacetimes as well as a tidally perturbed Schwarzschild solution 
(section 0. In these examples wc will see that the spacetimes that we would intuitively expect to be near equilibrium 
do indeed contain slowly evolving horizons. However we will also see that some fairly extreme situations also turn out 
to be slowly evolving. For example, wc will find that even the formation of a black hole during gravitational collapse 
can be a slowly evolving process. 

II. PRELIMINARIES 

We begin with a review of the definition of a slowly evolving horizon (section [ll All and a brief consideration of the 
units that we will use in our examples (section III Bp . 

A. Slowly evolving horizons 

A few preliminary definitions are needed before we consider that of a slowly evolving horizon proper. First, a 
marginally trapped surface (MTS) is a closed two-surface with flj-^) = and 6'(„) < Here 9(^£-^ = q°'^V a^-b and 
^(ra) = 9°''Va'T.& are respectively the expansions of the outward {€) and inward (n) forward-in-timc pointing null normal 
vectors to the two-surface and qab = gab + iant + Uait is the transverse two-metric on that surface. Throughout this 
paper we will cross-normalize the null vectors so that £ ■ n ~ —1 and tacitly assume that all MTSs are topologically 
S^. 

Then, a marginally trapped tube (MTT) is a three-surface of arbitrary signature that can be entirely foliated by 
MTSs 1^ and a future outer trapping horizon (FOTH) is an MTT which also satisfies 5n0{tj < 0. Here (5„ refers to a 
deformation generated by n and so the condition says that each MTS can be deformed into a surface with 9(^f^) < 
by an arbitrarily small evolution "inwards" . For a sufficiently small deformation 0(„) must also remain negative, and 
so this condition says that there are fully trapped surfaces "just inside" the MTT and untrappcd ones outside. Thus 
a FOTH can be regarded as a three-dimensional generalization of an apparent horizon ^. 

In defining a slowly evolving horizon, it is convenient to further restrict the scaling of the null vectors. To do this 
we label the foliating MTSs with a parameter v and choose the scaling and an evolution parameter C so that 

ya ^ ga _ ^^a f^^^ 

is tangent to the MTT and 

Cvv^l- (2) 

Thus V evolves the MTSs into each other and in particular the C characterizes how the area element y/q on the 
two-surfaces changes with increasing v. We have 

^vV^ -C0(n)\/1 (3) 

and so the area increases if C > (the horizon is spacelike), stays the same if C = (the horizon is null), or decreases 
if C < (the horizon is timelike). 

The value of the expansion parameter is directly determined by the behaviour of the matter and gravitational fields 
at the horizon surface. It was first shown in that if the null energy condition holds, a FOTH expands if the shear 

a^fj = (jaQb^c^d and/or the matter fiux Tab£°'i'' are non-zero. That is, it is spacclike and expanding if there is a flux 
of matter or gravitational energy across the horizon. Otherwise, it is null and so is an isolated horizon. This may be 
thought of as the second law for FOTHs and in fact it has recently been sharpened : if either a^^ 7^ or Tab£°'(.^ 7^ 
anywhere on a particular MTS of a FOTH, then one can apply a maximal principle to the partial differential equation 



This intuitive picture must be qualified. It is well-known that apparent horizons are slicing dependent. Similarly FOTHs are not 
uniquely defined and may themselves be smoothly deformed. Thus, while it is certainly true that a given MTT separates an associated 
set of trapped from untrapped surfaces, it will also be true that other equally valid MTTs and even f ull trapped surf aces will intersect 
it (though in restricted ways). These issues are still under active investigation. See, for example, Ia.l7i ri0lllllll2tll3l . 
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for C defined by 5\>9^£) = to show that C > everywhere on that two-surface. By the same reasoning, if C = 
anywhere on an MTS, then it is zero everywhere. Thus, there are no "partially isolated" MTSs and each MTS is 
either expanding everywhere or expanding nowhere [a, UM, UM ■ One consequence of this is that a FOTH can be cleanly 
split into isolated and non-isolated regions and the boundaries between these regions are two-surfaces of the foliation. 
We finally come to our definition |6|, |1J| : 

Definition: A region of a future outer trapping horizon H with v G ['yi,i'2] is a slowly evolving horizon 
(SEH) of order e if 



51. e= ^Cel^^Rn^l , 

52. the null vectors are scaled so that |V| = \/2C ^ e, 

53. \CvCoa\ < e/R% and l^v^^)! < e/Rjj and 

54. n, pp, |cr(")|2, and Tabn^n^ - 1/R]j or smaller. 

In the above Rh{v) is the arcal radius of the corresponding MTS, TZ is the two-dimensional Ricci scalar on those 
surfaces, uja = ~'7a"'cVbi?'^ is the well-known angular momctum one-form (or equivalcntly the connection on the normal 
bundle), and cr^^'' is the shear associated with the ingoing null direction n. 

Very briefly, these conditions may be interpreted in the following way. First note that even though both C and 
6'(„) explicitly depend on the scaling of the null vectors, the combination CSfn) docs not. For an isolated horizon it 

vanishes, while if H is spacelike we have C -pV? = — \fCj2Q (j^-^^fq^ where V = V /^/2C. Thus = CO'^^^^R^ is roughly 
the square of a normalized rate of expansion of the area element and clearly independent of the foliation parameter v. 
Condition ^is the key defining condition for slowly evolving horizons and in the upcoming examples we will often 
refer to e as the slowly evolving parameter. 

The second condition is not a geometric invariant of the surface. Instead it is a matter of picking a convenient 
labelling for the surfaces (equivalcntly a scaling of the null vectors) to reflect their slowly evolving nature. One way to 
ensure that this condition holds is to scale the null vectors so that 0(„) ^ 1/Rh- In particular for all of our examples 
we will choose = 2/ Rh (which is the value that it takes for Schwarzschild spacetimes with a standard scaling of 
n). 

One of the main reasons for choosing the foliation parameter in this way is to make the last two conditions easier 
to state. Then, condition ^says that the angular momentum one- form and 0(„) also change slowly with respect to 
such a parameter, while condition ^requires that conditions in a neighbourhood of the horizon not be too extreme. 

On applying these conditions to the evolution and embedding equations for FOTHs, it can be shown that an SEH 
also obeys dynamic zeroth and first laws of black hole mechanics. If one defines the surface gravity of the horizon by 

KV = -VntS/at (4) 

in analogy with the definition for event. Killing, or isolated horizons then 

Kv = K^"^ +0(e2) (5) 

for some constant k^^^ and function k^^\ Over any isolated segment of the horizon e = and one recovers the usual 
zeroth law. However, even when the horizon isn't strictly isolated, we see that the surface gravity only varies by order e 
across any MTS and, in fact, between MTSs during evolutions. This is what one would expect for a quasi-equilibrium 
state. 

The restrictions imposed by the slowly evolving conditions are not strong enough to fix a unique value for k^^\ 
However, at least for spherical symmetry, it is fairly easy to see that a natural scaling choice recovers the usual 
Schwarzschild value. This follows from the embedding equations which tell us that along a spherically symmetric 
MTT Q : 

C v0(„) = -nvOin) - + SnTabrn'' + C {SnTabn'^n'' + 9f^^/2^ . (6) 



Then if we scale the null vectors so that 6'(„) = ~2/Rh, we have 

1 

2R^ 



Kv = -^{1 + 8TTRj,Tab{r + Cn')n'') . (7) 
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Thus, if Tab = (as is the case for Schwarzschild spacetime) we find that k = \/2Rh- However even if it doesn't 
vanish completely, it is clear that we will have = 1/2Rh if Tabi^" + C'n!')n^ is small relative to the size of the 
MTS. In our examples we will see that this is the case when the horizons are slowly evolving. 

That said, in many cases such a "correct" scaling will not be obvious and then it is important to keep in mind that 
the SEH will still obey a zeroth law no matter what the choice. Further it can be shown that to second order in e it 
will also obey a first law : 

l«(o)^«,/ d2:,y~(T^,«;>+ 1 |^W|2\ (8) 



SttG dv Jh^_ V ^ 

where is a foliating MTS and aniv) is the area of that surface. This result recalls the standard physical process 
versions of the first law |l5l | , the first law for dynamic event horizons [l^ and the flux laws for dynamical horizons |3| ■ 

In cases where these horizons have an (approximate) axis of symmetry, it also makes sense to define their angular 
momentum and this can be incorporated into the first law in much the way one would expect. However as our examples 
will all be cither exactly or pcrturbativcly spherically symmetric, we will not consider the angular momentum terms 
in much detail here. 

A more detailed discussion of all of these matters may be found in . 



B. Units 



To build an intuition as to when a black hole will be slowly evolving we will often need to translate our examples into 
physical units. Since we are primarily interested in astrophysical processes, we take our standard unit of measurement 
to be one solar mass : 

Mq = 1.9 X lO^^kg. (9) 

Then, using the gravitational constant G = 6.7 x 10~^^ m^ kg~^ s~^ and the speed of light c = 3.0 x 10^ m s~^ distances 
and times that come out of our equations will be measured in units 

Rq = GMq/c^ w 1.4 X 10^ m and (10) 
To = i?o/c w 4.7 X IQ-^s (11) 

respectively. 

In some of the examples we will consider dust and have expressions for its density. Then, it will be useful to have 
a few density reference points for comparison. Note that 

Pq « 1.3 X I0\g/m^ w 1.9 X 10~^^/Ml , (12) 
p^^ w lO^kg/m^ w 1.4 X W^^/Ml and (13) 
pns ~ lO^^kg/m^ « 1.4 x 10-V^4 (14) 

are the mean densities of our sun, a typical white dwarf, a typical neutron star respectively. 

Finally for simplicity and consistency wc will typically study black holes that cither start or finish with solar mass. 
Such black holes are not fully realistic since astrophysical black holes forming from stellar collapse have masses of 
4 — 15A/0. However the conversion of our examples to more realistic masses is easy. For example to retool our 
examples with 6Mq as the standard mass, simply multiply the given times and lengths by 6 and divide any densities 
by 36. 



III. SLOWLY EVOLVING VAIDYA HORIZONS 



A. General considerations 

The Vaidya solution is the simplest example of a dynamical black hole spacetime j^- It models the collapse of 
null dust and is described by the metric 

ds^ = -fl-^I!!p.\dv^ + 2dvdR + R^dn\ (15) 
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where v is an advanced time coordinate, m(v) measures the mass of the black hole on a hypersurface of constant v 
and the infalling null dust has stress-energy tensor 

dm/dv 

Tab = 47rj^2 [dv\a[dv\b. (16) 
We consider spherically symmetric MTTs. Scaling the null vectors as 



it is a quick calculation to show that 



and n" = [0,-1,0,0], (17) 



1 / 2m(v)\ , ^ 2 

- 1-^ and (18) 



and so there is a FOTH located at Rh ~ 2m{v). Then from C\i6(^) ~ (since 6'(f) = everywhere on an MTT), a 
straightforward calculation shows that 



and so the slowly evolving parameter is 



dv dv 



^ Cel^iR^r ^ 8^ . (20) 



Then, conditions ^and 321will both be satisfied if dm/dv is sufficiently small. Given the spherical symmetry, ^and 
^are also easily checked and seen to hold in this case. Thus for Vaidya we have a slowly evolving horizon whenever 
dm/dv is sufficiently small. 

Further Tabi^n'' = Tabn°-n^ = and so from 10 it is clear that 

exactly. Then around a given value Vo we can Taylor expand to find : 

«v = f 1 - ^^A^vf + 0(e*)) . (22) 

This is in agreement with (jSJ (with k\ = 0). The surface gravity is constant over a given MTS and further only 
changes slowly between MTSs. 

In this case the first law actually holds exactly for any value of dm/dv. Specifically, both sides of equation © 
evaluate exactly as ^-Kdm/ dv. As we shall soon see however, this is a special property of Vaidya. In most cases this 
form of the first law only holds for slowly evolving horizons. 



B. Piecewise linear example 



Let us now restrict to a specific example to get a physical feeling for when Vaidya is slowly evolving. Consider a 
piecewise linear mass function : 

r 1 v<Q 

m{v) = I l + av Q<v <Vo . (23) 
[ 1 + avo V > Vo 

Thus, a black hole is irradiated with null dust for a finite period of time. For definiteness we assume the solar unit 
system discussed above and so this corresponds to a solar mass black hole that evolves to a (1 + aVo) solar mass black 
hole. 

Since there is no natural way to measure time along a spacelike horizon, a common problem in all of these examples 
is to decide on a reasonable reference frame in which to judge whether we would expect the horizon to be slowly 
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evolving . In this case, we make use of a fleet of observers who use rockets to hold themselves at constant areal radius 
Ro- Then, these observers will see a total mass of aVo fall past them during time 



T . £° ^1 - '-^dv < av. , (24) 

with the inequality coming close to saturation for very large Ro- 

Now, though these observers can't actually see the evolving horizon, let us consider the situation for a = 1/80000 
and Wo « 2.1 X IO^Mq (one second in standard units). Then, our observers would see aVo ~ 2.6 solar masses of material 
sweep past them in less then a second. However, at the horizon a quick calculation shows that 6^ = 80; = 10~'^. Thus, 
even in this fairly dramatic situation the horizon would be slowly evolving to order e « 0.01. 

IV. SLOWLY EVOLVING TOLMAN-BONDI HORIZONS 

A. General considerations 

Our next series of examples are Tolman-Bondi spacetimes and so consist of spherically symmetric clouds of timelike 
dust collapsing to form (or falling into) a black hole. These may be described in the following way (a more complete 
discussion of these spacetimes and their properties with particular emphasis on marginally trapped tube evolutions 
may be found in 0). 

We will consider gravitationally bound Tolman-Bondi spacetimes which contain an instant of time symmetry. That 
is, they contain a spacelike three-surface So with intrinsic metric 

ds" = ^\ ^, + r^idO'' + sin^ ed<P^) (25) 

1 — 2m\r)/r 

and extrinsic curvature Kab = 0. m{r) is a mass function defined by 

TO(r) = / dr{4:TTr'^)po{r) , (26) 

J o 

where po specifies the initial distribution of (stationary) dust on Sq. Thus on this initial surface the dust has stress- 
energy tensor 

Tab = PoUaUb , (27) 

where is the unit timelikc normal to the surface. Note too that the spacelike assumption means that r > 2m(r) 
everywhere and there arc no black holes in this initial data. 

Evolving this data forward in time, the dust falls inwards and one obtains the full four-dimensional metric 

ds' ^-dr' + ^^^^P^dr'+R\T,r)dn\ (28) 
1 — 2m[r)/r 

where t records proper time as measured by observers comoving with the dust. R{t, r) is the areal radius at time r 
of the spherical set of observers who had initial areal radius r and must be a solution of 



R'{T,r) = dR{T,r)/dr and the stress-energy tensor keeps the form Tab = p{T,r)\'aTVbT. 
In fact, this equation has an exact parametric solution given by: 



i?(r/,r) = rcos^ (^) and (30) 
T{v,r) = (?7 + sin?7) (31) 



.2. 
„3 \ 1/2 



\8m{r) 



7 



for a parameter 77 € [0, tt]. In our examples we will make use of this solution to generate the full spacetime evolutions. 
However, as an immediate application we note that it immediately follows that shells with initial areal radius r = 
will collapse to zero area at 



For all of the examples which we consider, this will be a strictly increasing function. In cases where it isn't, shell- 
crossing singularities occur (see for example 

For our analysis we need a pair of cross-normalized null vectors which we pick to have the form : 



ia = -v/1 - 2m{r)/r[dT]a + R'iT,r)[dr]a and (33) 
1 R'(t r) 



2^1-2m{r)/r' ' 2(1 - 2m(r)/r 
Then, it is a straightforward calculation to show that 



Rir,r) \^l-2m{r)/r J ^ ' 

and 

Combining the expression for with H29() it is easy to see that the trapping horizons may be found whenever 
R(t, r) = Rh = 2m(r) and on those horizons 

^(^)\h = -^<^- (36) 

Thus the surface i?(r, r) = 2m{r) is a FOTH. 

With this choice of scaling for the null normals we find that 

/ 2m(r)\ m'(r) , , 

C ^ 2 1 ^ — / ^ , 37 

V r ) R'iT,r) - m' r) ^ ' 



and so the first two slowly evolving conditions (^and will be met if 

EE Cel.Rl = (38) 



is sufficiently small. 

We have not given an explicit form for the parameter v. However, note that 

^^" = ^-2(l-2m(r)/r) ^^^^ 

and so when r 3> 2m{r) (which it usually will be in our examples) and the horizon is slowly evolving, C^t 1. The 
final two conditions (33 and are also easily shown to hold if e is sufficiently small. 
Next, consider the surface gravity. From {Tj) and (|27|1 we have 

^v-77^fl+ .,, f , w ■ (40) 



(41) 



2Rh V 2(1 - 2m(r)/r) 

Thus, when r ^ 2m(r), Ky will be 1/2Rh to order and it is also clear that k\> « ~C / {2R\) is small. 
Then, if as usual we assume that r ^ 2m(r) we have 

dm Ky rfflff 

dv SttG dv 

and 



V^r.''«' = 7^^77T4r^TTT-?^> (42) 



2 + C(l-2m(r)/r) dv 



and so the first law holds to second order as expected. 
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FIG. 1: The evolution of an initially solar mass black hole as dust shells of various masses collapse onto it. The bottom part 
of this graph shows the corresponding values of the slowly evolving parameter. It is clear that while all of the MTTs start and 
finish as slowly evolving, for larger mass shells they are not always slowly evolving throughout their evolution. 



Accretion 



For our first set of explicit ex amp les, we study accretion of a (spherical) dust shell of timelike dust onto an existing 
(and initially isolated) horizon |l8l |. To model this we start with both a Schwarzschild solution of mass nio and the 
Tolman-Bondi solution which evolves from a dust density distribution of 

on the initial surface Sq. Here, ^ is the total shell mass^, a locates its centre in units of Tq, and Tq characterizes its 
thickness : 86.5% of the mass lies within ±ro of aro, and 99.4% lies within ±2ro. 

Our desired spacetime is constructed from these two with the help of a bit of spacetime surgery. We remove the 
region inside r = 2mo in the Tolman-Bondi spacetime and replace it with the corresponding part of Schwarzschild. 
This can be done smoothly and gives a spacetime in which the dust shell collapses onto the black hole. The mass 
function of So is then: 

m{r) = rrio + 4tt df{?'^Po{f)} (44) 

J2Ma 

For definiteness fix the parameters of this spacetime as follows. Working in solar mass units, take mo = 1 so 
that the initial black hole has the same mass as our sun. Next, fix a = 4 and = 500 so that the dust shell has 
peak density at aro = 2OOOi?0 sa 2800km and almost all of the dust (99.4%) hes between 2000i?o « 1400km and 
3OOOi?0 w 4200km. Further, with these choices of a and Tq, w 2 x 10^^^/i and so for ^ 1 peak density is about 
that of a white dwarf. 

Then, with the help of the parametric solution given in (|31(l and (|30|l we can study the full spacetime. Specifically 
for fi = 0.4, 1.0, 1.6, 2.2, and 2.73M0 the evolution of the horizon is shown in Figure The top half displays its 
mass as a function of r while the bottom part shows the corresponding values of the slowly evolving parameter. In 
interpreting this graph, keep in mind that since r is given in solar mass units, the time span shown corresponds to 
about one second. 



^ The total mass will actually be very slightly smaller than fj, due to the excision of the region inside r = 2mo however for this example 
the difference will be negligible. 
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FIG. 2: The collapse of dust clouds of various initial masses to form a black hole. The bottom part of this graph shows the 
corresponding values of the slowly evolving parameter. While in all of the cases the MTTs asymptote into SEHs, note that for 
fi = 0.25 the horizon is slowly evolving throughout its formation and evolution. 

These values of /x were deliberately chosen to be on the boundary between slowly and non-slowly evolving horizons. 
It is clear that each horizon will be slowly evolving at the beginning and end of the accretion, this will usually not be 
the case during the peak inflow of matter. Note however that for fj, = 0.4, e < 0.01 throughout its evolution and so 
this case could reasonably be considered as slowly evolving throughout. This is true despite the fact that observers 
comoving with the dust would experience the OAMq of dust falling into the black hole in less than one second. 



For our next example we consider what would seem to be an even more dramatic situation — the gravitational 
collapse of a dust cloud to form a new black hole. Again we work with a Tolman-Bondi solution this time choosing 



where ^ is the total mass of the dust cloud, and quantifies its width : in this case 42.8% of the mass lies within 
To of the centre, 95.4% within 2ro and 99.96% within 3ro. 

To be concrete we pick vq = IOOOMq « 1400km and so this is a (spatially) small dust cloud. Further, with this 
choice the central density of the initially stationary dust is p(0) « 1.8 x 10~^°/x and so is somewhat higher than the 
average density of a white dwarf. We exphcitly consider the cases = 0.25, 0.5, 1.0, 2.0, and 4.0Mo. 

Then, the evolution as a function of r is shown in Figure [3 From the top half of that figure we see that horizons 
form between r « 20000 and r w 85000 (that is from 0.1 to 0.4 seconds after the collapse begins) and in all cases 
everything is essentially over by r = 175000 (before a full second is over). Turning to the bottom half of the graph 
we then see that while all of the horizons asymptote into the slowly evolving regime in the expected way, the higher 
mass solutions begin outside of this regime. Meanwhile, for yU = O.25M0, e < 0.01 throughout the entire collapse. 
Thus we see that under the correct circumstances horizon formation can actually be a slowly evolving process. 

Note too that again we have choosen our examples to be on the borderline between slowly and non-slowly evolving 
horizons. Lower mass and/or more dispersed clouds all fall even more firmly into the slowly evolving regime while 
more massive and/or more concentrated clouds will have larger values of e. 



Collapse 




(45) 



V. 



TIDAL EXAMPLES 



Finally, we (pcrturbativcly) break the spherical symmetry of the previous examples and consider the evolution of a 
black hole sitting in a slowly changing external field. Physically, such a situation would occur as a black hole of mass 
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M moves through an "external" gravitational field whose characteristic radius of curvature TZ as well as characteristic 
rates of change in time T and space C are all much larger than M . In certain regimes, one can split the total 
gravitational field into an external field and that part arising from the black hole and its motion. In particular, close 
to the hole one can describe the total field as that of the hole plus a perturbation from the external field while far from 
the hole the roles reverse. A specific example of such a configuration of fields would be a small black hole moving in 
a slow orbit around a much larger one. The formalism necessary to describe this situation was first developed over 20 
years ago [23, though we shall most closely follow the recent work of Poisson |23, 113 which pushes this formalism 
to the higher order that is necessary for our calculations. 



A. The Set-up 



In a little more detail, we begin with the external spacetime and the timelike geodesic 7 that the black hole will 
ultimately follow. Then in a comoving coordinate system the metric in some neighbourhood of the geodesic may be 
written as a perturbation of Minkowski space. In particular, at the origin of that system, coordinate time will coincide 
with proper time along the curve. The coordinate system will be based on the timelike unit vector u° up the geodesic 
plus a orthonormal spacelike triad e°, i £ [1, 2, 3], which is parallel propagated along 7. 

We assume that the external spacetime is vacuum and so the spacetime curvature is described entirely by the Weyl 
tensor Cabcd- With the help of the completely anti-symmetric Levi-Cevita tensor, this may be decomposed into two 
symmetric, trace-free, three-dimensional tensors: 

= ete''^CacMn''u'' and B,, = ^e^e'^e./'Cdebfu'uf , (46) 

which are respectively referred to as the "electric" and "magnetic" Weyl tensors. The perturbation expansion of the 
metric is then constructed from these tensors and their derivatives. 

To do this, we begin by making our earlier assumption about the magnitudes of quantities precise in the following 
way. First, the assumption about the radius of curvature of the spacetime says that the components of the electric and 
magnetic Weyl tensors arc on the order of (or smaller than) Tl~^. Next, with respect to geodesic time we have both £ij 
and Bij < TZ^^T^^. Similarly the homogeneity assumption says that on extending the tetrad into a neighbourhood of 
7, spatial derivatives of these quantities such as C ^i^ij ^ TZ^^C^^. For notational convenience we reduce these three 
parameters to one, defining X = min{7?., £,T} and consider the expansion in powers of X"^. Thus, the components 
of the electric and magnetic tensors are < while the derivatives of these quantities are < 

Now, in Fermi normal coordinates based on the freely falling {u°',ef} tetrad the metric may be expanded in the 
following way. First, since the tetrad is orthonormal, the metric is Minkowski to zeroth order. Next at dipole {X~^) 
order there are no corrections since the frame is freely falling. The first corrections occur at quadropole {X~^) order 
where they are are constructed from the £ij and Bij. Terms at octopole {X~^) order depend on the temporal (along 
7) and spatial (in the ef directions) derivatives of these quantities. Higher order terms similarly depend on higher 
order derivatives though they will not concern us in this paper. 

If we insert a mass M X Schwarzschild black hole into the spacetime and set it to travel along the worldline, we 
can characterize the induced perturbations near the hole {2M < r ^ A") in terms of those same quantities. Specifically 
in ingoing Eddington-Finkelstein type coordinates [23l |. 

ds^ = " " T" ~ '^^™) " '^^^^^ + 2SgyAdvde^ + {r^flAB + Sqab) d0^d0^ , (47) 

where JIab = diag[l, sin^ 9], and A,B & {1, 2} with 9^ — 9 and 6*^ = (j). The tidally induced perturbations then take 
the form: 

3 

5g™ = r^eif ^ - y (ea^^ - 63^°) + 0{X-^) (48) 



S9AB = -V(e7fli3-&76^B)-f^M^B-^8e^B)-^MB-&9SAB) + 0(A'-4), 



where the exact definitions of the various £ and B terms are given in appendix The q and o superscripts indicate 
the multipole classification (and hence expansion order) of the various terms while the over-dots indicate derivatives 
with respect to v (and so lower the order of dotted terms by a factor of X). The full expressions for the e„ and 6„ 
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are given in appendix though here we note that if M = (ie. no hole) the e„ and &„ are all unity while if M 7^ 
they are functions of r alone. These functions go to unity as r — > cx) and vanish at r = 2M except for e7(2Af) = 1/2, 
eg{2M) = 1/10, 67(2M) = -1/2 and bg{2M) = -1/10 (though the derivatives arc, in general, non-zero). 



B. Locating a horizon 

In general; locating an MTS is not a trivial task and indeed after a glance at the metric (|47|) . one might suspect 
that that will be the case here. In fact, such worries are unfounded. It turns out that to the order in which we are 
working, r = 2M foliated by u = constant surfaces is a FOTH. 

The see this, first note that the induced metric on a r = 2M and v ~ constant is 



(49) 



qAB = UP ilAB - i:M{£1b + ^Ib) ^ T^HSab + B\b) + 0{X-') 



and so the corresponding area element is 



(7 = 4Af^ sin6' + 0(A'"*), 



(50) 



which follows from the fact that the perturbations are tracefree with respect to the unperturbed (spherical) metric. 
Next, up to rescalings, the future outward pointing null normal vector field to these surfaces takes the simple form 



[dr]a ^ r 



d_ 

dv 



thanks to the vanishing of the majority of the e„ and 6„ at r = 2M . Then 



1 



Hi) 



1 d 
= + 0(^-4) 



(51) 



(52) 



by direct calculation. 

Shifting our attention to the inward null expansion, the future inward pointing null normal vector field to the 
two-surfaces (cross-normalized so that £ ■ n = —1) is 



-[dv]a <^ n" 



d_ 

dr 



+ 0(^-4), 



and so 



Finally another direct calculation shows that 



SJ 



d ^ 

dr 



M 



+ 0{X~^). 



,)+0(A-4) 
+ 0{X-^). 



(53) 



(54) 



The in the last line indicates that the metric can be used to calculate this quantity to higher order (actually 
X~^) than shown, but this has not been done since all that we needed to show was that 5nd(() < 0. 

Thus, combining the results of these three calculations we see that the r = 2M surface foliated by u = constant 
surfaces is a FOTH as asserted. 



C. Horizon properties 

Next, it is straightforward to sec that this FOTH is actually a slowly evolving horizon. First from equations and 
(I50|l it is immediate that the expansion parameter C is at most of order X""^. Combining this with the expression 
(|53|l for 6'(„) it is easy to see that conditions ^and ^hold. 
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Checking the other slowly evolving conditions requires the calculation of a few more quantities. Since the spacetime 
is vacuum, the matter terms clearly vanish while the other quantities are: 

= (^Ib + + 8Af^ (-y^^B + + 0{X-') , (55) 

= -^M' {£\ - B\) - ^M^ (\Q£\ + 'd£\ - 36^) + 0{X-^) and (56) 

^^^2+^(^'')^ (57) 

where the "w" in the expression for the Ricci scalar again indicates that the metric actually defines this quantity to 
higher order than indicated. Then it is easy to see that conditions ^and ^also hold and so these tidally distorted 
black holes have a slowly evolving horizon of order e = A""^ at r = 2M . Further, though we will not show the explicit 
calculations here, it is not hard to show that the standard two-sphere rotational Killing vectors remain as approximate 
symmetries (as defined in 0) of the perturbed horizon. 
We can then consider the mechanics of this horizon. First 

«v = ^ (l + y ^^'^') + Ol-^"') • (58) 

Thus, the surface gravity is approximately constant over each Hy and further only changes slowly in time as would 
be expected for a quasi-equilibrium state. Note too that the expression for the surface gravity matches that found in 
|23| (though this isn't especially surprising for this spacetime since the event horizon overlaps this FOTH to the order 
in which we are working). 

The first law also holds. Since C is at most of order = X^^, the left hand side of equation JSJl vanishes to this 
order. To evaluate the right-hand side we need an expression for (t\^ which we calculate as 

^AB = VQAB = -^M' {sIb + B\b) - {^AB + B°ab) + ©(A^^) . (59) 

As in the expansion calculation H52|l . the Lie derivative is simply evaluated as the time derivative with respect to 
V. Then given the form of H49|l and the assumption that taking derivatives with respect to v lowers the order of 
metric quantities by a factor of X, the above expression follows directly.^ Thus, the right-hand side of that equation 
evaluates as: 

^ d^x^qai'^a^'^^" = (4^^ + B,,B^^) + ^A/^ (^.r^^ + B.,,B^^'^) + 0{X-^) , (60) 

where the computation makes use of equations H5U|1. (|59|1 . and (jA5IAll|l and is straightforward but quite lengthy. 
Since both of the left and right hand sides of |(HJ) are zero to order X^^ the first law holds. 

More interesting though is the result (|H5|l . As it is the "square" of a quantity known accurately to order X~'^ we 
can trust it to order X~^ . Thus, we see that in this case the first law allows us to calculate the rate of change of the 
mass even though we cannot see those changes directly in the approximation that we are using. Note too that this 
result matches those found in [2^ |2^ though with the difference that our calculation gives a snapshot rather than 
time averaged value for dm/dv. 

Finally we again consider the physical situation under which this is a slowly evolving horizon. Taking M = A/q the 
time scale of the changing fields must be significantly larger than Tq. Taking T = IOOTq and assuming that A" ~ T 
this gives ^ ~ 10-12 ^ lO-^A/g/sec. 



VI. CONCLUSION 



This paper provides concrete examples of slowly evolving horizons and demonstrates that they correspond to 
situations that one would intuitively consider to be slowly evolving. Thus we saw that Vaidya and Tolman-Bondi 



While we have only claimed fourth order accuracy in 1511 , in fact horizon-locking coordinate choices have been made the mean that this 
shear may quantity may be calculated directly from the m etric to thi s higher order. This is confirmed by a gauge-invariant calculation 
based on the even-parity and odd-parity master functions l23l l24t l25l . 
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black holes with sufficiently small mass flows are slowly evolving as is the tidally perturbed Schwarzschild black hole. 
In the Tolman-Bondi examples we studied horizons that start out isolated, transition through the slowly evolving 
regime as an infalling dust cloud begins to accrete, pass out of it as the bulk of the matter hits, and then return to 
that regime as the last dregs of dust fall in and the horizon asymptotes towards isolation. Similarly, in the aftermath 
of a black hole formation from a dust cloud collapse, the MTT is slowly evolving as it approaches isolation. A bit more 
suprisingiy, we saw that in some cases the formation of a black hole during gravitational collapse can be characterized 
by an SEH throughout its entire evolution. 

It also became clear that, at least in the spherically symmetric case, we need to pay careful attention to mass scales 
in intuitive judgements as to whether a particular horizon is slowly evolving or not. For approximately solar mass 
black holes accretion has to be truly spectacular in order to have a non-SEH MTT. By contrast, for a supermassive 
IO^A'/q hole the characteristic time scale would be over an hour and the characteristic density would be ~ lOOOkg/m'^. 
Thus, in contrast to smaller black holes where matter of density somewhere between that of a white dwarf and a 
neutron star is needed to move an MTT out of the slowly evolving regime, here only densities about that of liquid 
water are needed. 

Given the orders of magnitude involved in these examples, it seems like that the intuition that we have gained here 
may apply, at least qualitatively, even away from spherical symmetry. Thus, the astrophysics of stellar mass black 
holes is probably largely "slow" (apart from especially dramatic events like black hole mergers, neutron star/white 
dwarf captures, or some formations) while that of super-massive black holes SEHs will often not be in that regime. 
For example, one suspects that the capture of even a small star by a super-massive black hole might break the SEH 
conditions (and it might even cause a horizon "jump" as discussed in |3,[la|)- We hope to see this conjecture tested 
by future numerical simulations of more realistic astrophysical situations. 

In developing any theory it is very useful to have specific, easily handled, examples in mind against which one can 
test one's ideas. This paper provides some of these examples for slowly evolving horizons and so helps to form both 
our mathematical and physical intuition about the quasi-equilibrium states of black holes. 
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APPENDIX A: IRREDUCIBLE TIDAL FIELDS 

To octopole order, the terms appearing in the metric are |23| |: 







cq 




cq 


= "^PAPB^ab + QAbS'^ , 


B\ 


= p^e.jk^^B'^in^ 


B\b 


= ''-I^(APBf^3f'B^^ 


£° 




CO 

^A 




CO 

^AB 


= 2p\p>g£,,k^*' + qAB£\ 


B\ 


= lp\e,jk^'B''i^n'n"\ and 


B'Xb 





(Al) 



where 

17* = [sin 9 cos 4>, sin 9 sin cf), cos 6] 
is the standard radial vector in Euclidean IR'^ and 

£i]k = {D^£jk + Dj£kt + Dk£ij) and 

B^Jk = \{D,Bjk + DjBk^ + DkB,j). 



(A2) 

(A3) 
(A4) 
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Each of these terms may be expanded in terms of the appropriate spherical harmonics, but for our purposes we 
note that direct integrations over the unit sphere with dfl ~ sui9d9d(f> give: 

j^ne^^^'Els = (^l^^£,^r\ (A5) 

dnB'"'''B\s = (^7t)b,,B'\ (A6) 



S2 



5 

dnr^^'E^^B = (§^)f.J;.r^^and (A7) 



jjnB"^^B\s = (^-^n^B,,,B^^\ (A8) 



AU other combinations are zero 



dn 8'^'^^B\g = / dn £°^^B°AB = , (A9) 
I dn £'*^^£°AB ^ I dn B'^'^^B^B = , and (AlO) 
dn8'^^^B°AB^ I dn£°^^B\B = 0. (All) 



APPENDIX B: RADIAL METRIC FUNCTIONS 



The radial functions appearing in the tidal metric components ()48ll have the form 



ei = P 

62 = /[I + 4^(5 + 121nx)- 31^(9 + 41na;) + 3|, 

64 = / 
65 



3 

42? 



/ [1 + ;f (13 + 121na;) - A - A - A 

3x) 



66 = /(l-l 

67 = 

68 = l + i(4 + 31n:r)-^-J,(7 + 31nx) + 3f, (Bl) 

69 = / + To^^ 

bi = f 

be = /(l-^) 

bs = l + ^(5 + 61na;)- A- A(2 + 31na;) + -^ 

bg ^ f ~ , 

where x = r/{2M) and / = 1 — 2M/r. Their derivation is non-trivial and there is a (coordinate) gauge freedom in 
their final form. The reader is referred to [2^ for details. Here we simply note that the freedom has been put to good 
use to ensure that many of these functions vanish at r = 2M. In fact the only non-zero ones are 67 = i, eg = j^^j, 
&7 = — i and 69 = —j^- This greatly simplifies our calculations. 
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